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ABSTRACT 


A decrease of one to six gauss (depending on magnetic 
field orientation) in the coupling constants of the CH.CO, | 
radical formed in X-1irradiated zinc acetate dihydrate has 
been observed over a temverature range of about -60°C to 
+30°C. Calculations of dipolar and Fermi contact intéraeerem 
based on a model of internal rotation of the methvlene grounv 
about the C-C bond have shown a small coupling constant de- 
crease on the order of 0.19 gauss over a temperature range 
of -150°C to +90°C. The major effect has been shown to be 
due to a spin relaxation mechanism. The effect was calculated 


using Monte Carlo techniques, the results of which were con- 


firmed by experimental data. 
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i. IN@RODUCTION 


Study of the election spin resonance (esr) spectrum of 
Single crystals damage by radiation have led researchers to 
the identification and structure of a number of paramaagnetic 
eeeeres (Refs. 8, ll, 15, 13, 21, and 22]. Radicals may be 
generated by irradiation of single crystals at nearly any 
temperature where the crystal remains intact. However, 
irradiation at one temperature will not necessarily produce 
the same radical or combination of radicals as that at an- 
Other temperature. Hence, by selection of the anvropvriate 
temperature, it iS possible to selectively study several 
radicals which may be generated in the same system. 

A considerably lesser amount of research has been done 
concerning the temperature dependence of the observed esr 
Spectra for given paramagnetic species. Crawford [Ref. 4], 
in his study of the zinc acetate dihydrate crvstal and the 
paramagnetic species, CH.CO., 


Jom 


temperature effect in the sum of the coupling constants. 


formed therein, noted a marked 


The effect ranged from a one to six gauss change (for an av- 
erage coupling constant of about 45 gauss) over a temvnerature 
range of 298°K to 133°K depending on the external magnetic 
Field orientation. Moss [Ref. 14] calculated the temverature 
effect in the esr spectrum of methyl radicals by attributing 
the effect to out-of-plane vibrations. However, he was unable 
to account for the magnitude of temperature coefficient bv 


applying his method to the benzene anion. In a study of the 


a 


hydrazine positive ion, NoHy Falle [Ref. 7] observed an in- 
crease in the ratio of nitrogen hyperfine splitting to nroton 
hyperfine splitting, Ja/ay | with increase in temperature. 
The effect was attributed primarily to a decrease in the vro- 
ton hyperfine splitting, a, as caused by torsional (rota- 
tional) displacement of an amino group about the N-N bond. 
This torsional motion causes an introduction of positive spin 
density at the protons leading to an increase in proton hy- 
perfine coupling, ane with increased torsional motion as ex- 
perienced with increasing temperature. Schrader [(Ref. 17] 
has theoretically (based on the procedure of Schrader and 
Karplus [Ref. 18]) predicted the counling constants for the 
methyl radical over a wide range of temperature from zero to 
500 degrees Kelvin, and has assigned temperature coefficients 
thusly which are in good agreement with Moss' results [Ref. 
14]. Zlochower, Miller, and Fraenkel [Ref. 23] also applied 
the theory of Schrader and Karplus [Ref. 18] to the methvl 
radical and observed experimentally a dependence of the pro- 
ton hyperfine splitting with temperature on the order of 
0.002 gauss/degree in a 50°C range extending above and below 
room temperature. Work done on the isotropic proton hyverfine 
interaction in the esr study of t-electron free radicals has 
shown that out-of-plane vibrations should have a contribution 
to the splitting as a result of direct overlap between the 
proton and the unpaired electron in the tm orbital [Refs. 13 


and 14]. This contribution would be expected to increase with 


a2 


increasSinc temperature leading to a corresponding decrease in 
the coupling constants with increasing temperature. 

Crawford [Ref. 4] gave no detailed expolanation of the 
marked shift in coupling constants observed for the CH,CO. 
radical of zine acetate. It becomes the dual vurpose of this 
work to verify the effect observed and then to develop a model 


by which the phenomena may be exvlained. 


I 


Il. BACKGROUND OF THE CH.CO, RApECAL 


Tolles, Crawford, and Valenti [Ref. 21] have recently 
published work on the epr (esr) study of zinc acetate dihy- 
drate single crystals. Almost simultaneously, Ohigashi and 
Kurita [Ref. 15] published a similar work. 
7C0> radijeam 


to exist in X-irradiated single crystals of zinc acetate. 


Both groups of authors have identified the CH 


The radical is stable for hours at room temperature with a 
half life of about ten minutes at 80°C. Figure 1 shows typ- 
ical esr spectra generated by the irradiated zinc acetate 
system. Coalescence of the typical four absorption dine 
spectrum observed at low temperatures into the typical 1l:2:l 
three line spectrum at high temperatures as indicated in Fig- 
ure 2 is evidence of internal rotation of the methylene group, 
CH., about the C-C bond of the radical. This phenomena and 
determination of the barrier to internal rotation is discus- 
sed by Tolles et. al. [Ref. 21] and Ohigashi and Kurita 
(Re fees |e. 

The crystal structure of zinc acetate dihydrate and the 
a, Db, Cc cell axis system as related to am a*,b, ¢c (x Fuga 
orthogonal axis system within the crystal is shown in Figure 
3. The relationship between the two axis systems 1s given bv 


the transformations 


14 


High Temperaturt 
[Oe CE rum 


spectrum 


| 
| 
Low Temverature | 
| 
| 


Figure 1. Characteristic esr snectra at high and 
at low temperature. 
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Lows Lenpeiareui.e 
Sjolsren wig) Sic 


Coalescence as 
Temperature 
Increases 


High Temperature 
Spectrum 


Figure 2. Coalescence of absorption lines of 
esr spectra with increase of 
LeEMpeLracure . 
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Figure 3. 





The zine acetate dihydrate crystal and 
the orthogonal crystal axis system (xyz 
or a*bc) and its relationship to the 

crystalline axis system (abc) {Ref. 4]. 


yy 


&@ = sinkx + cosBz 
b= y 
CZ 


where in Figure 3, x 1s also referred to as a*. 

The ein o> radical is shown in Figure 4. The orientation 
shown is planar, each C-H group having assigned to it an 
orthogonal axis system Xi Yas Zi or Xo, Yor Z. wherein the 
Ya and Y, axes are coincident and fall along the om orbital 
of the carbon (C.) atom. These directions correspond to vec- 
tors u, Vv, and w in the a*,@s, C Cayciaal axis System. 

The following parameters of the CH.CO, radical as pre- 
sented by Tolles et. al. [Ref. 21] and Ohigashi and Kurita 
[Ref. 15] are pertinent to studies in this paper. 

Principal values of the hyperfine coupling tensors” 

Hy tensor <-26.80 MHz 
=59.59 
-92.15 


H. tensor -26.41 


2 
ea) ye 2) S 


69 2 cms 


+principal values are designated to lie along the resnec- 
tive C-H, bond, the De orbital axis of the carbon atom in the 


CH, group, and in the CH, plane perpendicular to the respec- 


tive C-H bond as shown in Figure 4. 


I) 


O; 


Figure 4. 


Vek, 





LOCaAe1ONS Obeene Nypewrine, axcom) maene 
CH,CO, radical (planar configuration) 


[Ref. 4]. 


AAD 


Direction cosines with respect to the crystal axis 


éyetem- 
Hy O.547 Ou 22 05423 
0.768 -0.634 0.088 
=O), JOC 2/6 GeoC2 
Ha 0.048" °-0-08)1 “=-0799'6 


0-268 —-0-654 0.088 


0 2638 0. 769 "—-0 ster 


The barrier potential to internal rotation of the 
methylene group about the C-C bond 

5.8 kcal/mol 
The average rotational frequency factor 


Awe xX 1072 seem 


The direction cosines are expressed as row vectors. 
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III. PROBLEM FORMULATION 


To qualitatively understand the effect being studied, 
reference is made to Figure 1, page 15 which shows a typical 
first derivative esr spectrum of the CH.CO, radical at high 
temperature and at low temperature respectively. The esr ex- 
periment basically measures the interaction of the election 
magnetic moment with an applied magnetic field, H. Hynerfine 
structure is due to the free election interacting with the 
two methylene protons. The coupling constant, Arps isa 
property of the radical dependent on orientation in the mag- 
netic field and reflects the magnetic field of a given nu- 
cleus, Hi, experienced by the electron. 

Figure 5 shows the splitting of energy levels in a maq- 
netic field due to electron and nuclei interactions for the 
CH.,CO, radical. The four vertical lines represent allowed 
energy transitions giving rise to the basic four line spec- 
mmam Of Figure lL. 

When an oscillating field perpendicular to the external 
magnetic field H is applied, transitions are induced pro- 
vided the frequency, v, is such that the resonance condition 


hv = g8H 


is satisfied, where g is the spectroscopic splitting term and 


8B is the Bohr magneton. Hence it can be said 
AE. = Av. = gBH, 
or 
H. <AE. 
i ai 


onl 


+] 


t+ 
Il 


Z 
aot : 
j i 
s=+5 I5=+1 
i | =—] 
L ( 
! 
| 
2 Ja 
» § | 2 
a4 a 
bt ft ot I, 1 
rt 
= ( 
Ij=-l | | | 
I px 
7 1 4 T5=-1 
: T,=+1 
Electron Zeeman : 
Splitting . 
Hyperfine Splitting 
Due to Hy I5=+1 
Hyperfine Splitting 
Due to H. 


Figure 5. Splitting of Energy levels@in@. magnetiestield @ue togelectron 
and nuclei (Hy and H,) interaction in the CH,CO. radical. 


Ze 


Fach resonance (absorption) laine of the spectrum is a repre- 
sentation Of a transition enekay AE proportional to the 
magnetic L£1ela, Hoy causing resonance. Hence, the difference 
between two spectrai lines in gauss 1S a representation of 
tue diifrerence In transition energy which eserr 1S a measure 
of the hyperfine interaction and is called the couvling con- 
Stant. Experimentally the "distance" between absorption 
fines is directly measured in gauss to determine the relative 


coupling constants, ny and Ay 1 OY the sum of the coupling 
1 2 


constants, a, +a seen in Figure i. 
ip Z 


As previously mentioned, Crawford [Ref. 4] in his work on 
Zinc acetate dihydrate noted that at a given orientation of 
magnetic field, the esr spectrum showed a one to six gauss 
decrease in the sum of the coupling constants as the tempera- 
ture was increased from 133°K to room temperature (Figures l 
and 2). In section I studies of temvnerature dependence of 
radical spectra have been discussed. These studies have con- 
tributed similar temperature effects to out-of-plane vibra- 
tions. Studies by Horsfield [Ref. 12] and Heller [Ref. 10] 
of methyl groups have shown internal rotations of the qrouv 
to occur about a C-C bond. The temperature dependence of the 
eaeco., radical has been attributed to motion of the CH. qroup 
consisting of rapid reorientation around the C-C bond above 
room temperatures being quenched at low temverature [Refs. 15 


and 2] |= 


Ze 


It 1s the purpose of this work to develop a mathematical 
modei, consistent with exnerimental data, by which the tenner= 
ature effect in the CHAO; radical may be vredicted for anv 
given orientation of the zinc acétate crystal in a maanetic 
field. 

It will be seen in thediscussion of results’ that the Gb- 
served effect does not depend solely on one particular vnhénom- 
enon but rather 1s duewbo a combinationWof severa: vhememenda. 
Hemee the problem is formulated into toummeffects, rotamer 
Of the methylene group, hyperconjugation, flan or vibratiege 
perpendicular to the plane of the radical, and nuclear Same 


relaxation. 


Zoe ROTATION OF THE METHYLENE “GROUP 
1. The Restricted Rotor Model 
The CH, (methylene group) acting as an Gnternally 
restricted rotor representing the exchanging of the two wro- 
tons 1s the model used in explaining the temperature effect. 


The Hammeitonian. fOr a Beckricted fotor 1s 


(l=eos2. 05) 


NS 


assuming the radical to be planar and the angular momentum, 
L,, to be quantized along the Z axis (the C-C bond). In the 
Hamiltonian equation, B 1s a constant of the system, V 1s 
the activation energy or barrier to internal rotation, and 9' 


imac angie Ge rotation Of “the CH Mrotor Gabout the C-—C bond: 


2 
The Schroedinger wave equation for the restricted 


rotor may be solved by such standard vrocedures as given in 


24 


Eyring, Walter and Kimble [Ref. 6], page 358, to yield wave 
functions which in the limit of the free rotor become: 


oe ies Ae? 


| | iL 
where A is a normalization constant equal to —-. 


V2T 
Tae set Of basiS ELUNCLIONS are venoscen erence 


l —aainoe 


MWeere m runs from minus to plus infinity. 

Elgenvalues of the Hamiltonian yield the possible 
quantum mechanical energy levels of the system. These values 
are established by diagonalization of the Hamiltonian matrix, 
the diagonal values being the eigenvalues. 

The Hamiltonian matrix elements are determined 
through the following operations: 

On-diagonal matrix elements are represented as 


2 


sm! | BLS e+e | ms = <m' |H'|m> = (Bm™+s) 6 


2) mm ' 
where 6 = 0 unless m = m'. 
The diagonal elements are numerically evaluated from 
the knowledge of V, the barrier potential, and B, 
where 

Ve= oc ey met 


y= he 
817° 1 


h = Planck's comewant 


I = the moment of inertia of the 
rotating CH, group. 


oS 


For conven.ence of calculation, 


- 
8721 


where 
oO 
gor = B40 55307 10° MHz amu Ae 


Off-diagonal elements of the Hamiltonian matrix result 
from the perturbation term 


V 
-5 cosoé 
and are represented by 
<m' | 4 cos¢'|m> 


Evaluation of these integrais shows that the value is 
zero unless 

m! ="mr2 or m=2 
The non-zero off-diagonal elements are then expressed 


in termGwer themeascis funceren as 





27 _~1imo' imo ' 
-¥ ( =—-—— cos26' = dg’ 
0 V2T V2T7 
Using the identity 
: t ee t 
Gos 20) )= see +e Bag ) 
the integral is evaluated to yield 
V : 
Om wien © = M+Zs0r ™m' —sM-2 etne Oft-en agonal walues 
become 
_V 
4 


26 


The calculated matrix elements are placed into the 
secular determinant which is diagonalized to vield the eigen- 
values and eigenvectors of the system. Solution of the deter- 
minant was performed by IBM 360 comvuter. 

2. The Harmonic Oscillator Model 

The limits of a restricted rotor model are the free 
rotor at one extreme and the harmonic oscillator at the other. 
Considering the first several excited states above ground to 
be significantly populated at the temperatures used in this 
study, an approximately constant difference between levels of 
the restricted rotor was seen to exist (Section V-B). 

The nearly constant difference in eneraqy levels suaq- 
gests the approximation that the system can be represented by 
a harmonic oscillator model, the energy levels of which are 
represented by 

1 (nt+5) hy (A2-1) 
where n runs from zero to infinity. If this approximation 
were to prove valid, further system analysis would be made 
eaSler. 

The restricted rotor model yields eigenvalues in units 
of E/h, hence some rearrangement of equation A2-1 is necessary 
to compare the two models. A solution of the Schroedinger 
equation for the harmonic oscillator by procedures of Eisberg 


feet. 5] shows 


- 


Vx (A2-2) 


Equation A2-1 may then be rewritten 


E = (nei =) hy (> . (A2-3) 


Z| 





Direct substitution with proper units as used in section A l 
y1ireids 


E/h = 2(n+5) VI8. 159208 (A2~4) 


From which eigenvalues may be determined in units of terahertz 
(THz) consistent with the units of the restricted rotor model. 
(One THz equals Toe varie 
So. EVALUAEEON OF @=COlUpe Tm Ccometrents asa Function oF 

Temperature 

The harmonic oscillator model was used to develon the 
values of coupling constants expected to be seen experimentaliv 
at any given temperature due to rotation of the CH, Grouwe 

THe ELteccive Spin HaMmrecontan for the tnteractiom 
of an electron and two nuclear spins in a maaqnetic field for 


the case of quenched angular momentum 1s 


dq = Bii-g-8 = ero mns 1-H + 
The first term represents the electron Zeeman energy, the 
interaction of an unpaired electron with the external magneieme 
field H. 8 1s the Bohr Magneton, S is the electronic spin 
equal to “s and : is the spectroscopic splitting tensor. The 
second term 1s the summation of the nuclear Zeeman energies 
where g, is the nuclear g factor, B is the nuclear magneton, 


and I. 1s the spin of the aaa preton. the third term Yvermes 


Senltememe Chergy Of Nyperiinemmuiveract: on between electron 


7 


>> 


— 


Softies, ana Nuclear spin, i. 1, 1s the hynerfine inter- 


ecelOn Censor. Equation AS-i, be remertten 
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ay se 
H= @h-g-6 - gee) or oe (A3-2) 


where the parenthetical term is the cftective magnetic field 
seen by the i proton. Solution of equaticnes3-2 teadsmeo 
(among other parameters) the coupling constants of the spin 
System. 

When an electron and aypYoton are placed an the anelu- 
ence of a magnetic field, the spins of each tend to line up 
parallel to the field creating dipolar counling which is 
represented in the effective magnetic field term of equation 
A3-2. If the orientation of the magnetic field is changed 
relative to the proton and electron as caused by the rotation 
of the CH, grouv, different dipolar counling occurs, and the 
Hamiltonian solution yields different values of the counlina 
constants. Hence, the coupling constants ere a function of 
the CH. anglewoLisrotation, oO” . 

The Hamiltonian is a representation of the energy of 
the spin system. There will be a certain probability ae 
that a given energy level will be found at a given angle of 
CH. mocation, @'. In addition, there will be a Statistical 
distribution of electron spin population in the existing en- 
ergy levels as determined by Boltzman distribution factors. 

Based on the foregoing considerations the general re- 
lation to be solved to yield the average coupling constants 


for a given temperature is 
ie -( Agr Py Wy de a2) 


ae) 


where <A>,, is the coupling constant averaged over all orienta- 


bh 
tions for a given magnetic field orientation at a given tem- 


perature, Tea. A is the counling constant as a functlonieg 


db i] 
angle of rotation, $', of the methylene groun about the Gee 


bond axis. P ismthe pmebabiieity distribution acter. Wa 


b! 
is the Boltzman statistical weighting factor. 
a. Eealluiaerone o£ ae 

An existing comouter program, SSESR (Solid State 
Spectrum ESR), was used as a starting voint to calculate val- 
ues of Ages As originally written by Crawford [Ref. 4], 
SSESR computes the esr spectrum of a radical tranned in the 
solid state assuming the angle $' to be zero, i.e. the nlanar 
model. 

The theory behind the SSESR nrogram is based on 
the development of Box et. at [Ref. 2]. If the comvonents 
Cus S perpendicular to H are ignored in equation A3-1 (a qood 


approximation [Ref. 2]) the energy levels of the snin svstem 


can be written as 


~ 


E = 8gHS,, - UPnP Tiny (A3-4) 


where g is the observed g value, S,, is the electron spin com- 


H 


ponent along H, and Tey is the comnonent of nuclear spin along 


aL the vector along which the nuclear spin is quantized. The 


vector hy may be written 


h Ss S CU A +V.B B.+w.C ) 
i Hoe ee a 


; (A3-5) 





Jo ta 


In A3-5, A, B. and Cc. are princival values of the hyperfine 


-> 


; > > -> > : 
COuUpnENG ECNSOLG ,.24. The vectors Ul, v., and’ w., are unite 
a si at 
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vectors denoting the direction of the two sets of princinal 


axes corresponding to A., Bay and Cc. respectively, as seen in 


Figure 4. The terms Os 8, and y, are the direction cosines 
of the external field, H, with respect to these axes. For 
= ale 
te 
the CH.,CO. ieaGune acl Sus Toy and Toy are each ix. 


Incorporating all allowed combinations of spin 
states, solution of equation A3-4 yields 16 electron transi- 


tions in terms of 2AHgGBG/G,8 represented bv 


h h (A3-6) 


| = | | 
ant) DS, 
where the 16 transitions are seen by the te choices of sign 
and where Paya) and |hy ¢_y | are the magnitudes of the vec- 


tors h. corresponding to Sy = ri and -5 resnectivelv. The 


2 
computer program, SSESR, calculates the relative intensities 
and frequencies for the 16 transitions of the pvlanar radical. 

Basic input data include the number of nuclei in 
the molecule or radical, magnetic field maaqnitude, direction 
cosines relating the nuclei principal axes to the crvstal 
axis system, the principal values of the coupling tensors for 
each nucleus, the direction cosines of the magnetic field, 
and the nuclear Zeeman frequency. Output yielded transition 
frequencies, transition intensities, and transitions in gauss 
from the center of the spectrum, all assuming an isotrovic 
g tensor [Ref. 4] or a value of g = 2.0023, that for a free 
electron. 

Several modifications to the existing program 


were made to incorporate the harmonic oscillator model. The 


Sul 


complete program, PREDICT II, appears in Appendix A. All 
computations were performed in double precision. 

(1) Subroutine sBeDecs.. The first modi ficatver 
involved an additional subroutine to compute the three direc- 
tion cosines of the magnetic field for a given input of CG vam 
o, the angles through which the crystal would be rotated in 
the magnetic field. This subroutine was called HDCOS and was 


based on the relations 


Uy = -cos 6 
U5 = -sin 6 sin 6 
U. = iba & wees uw 
(2) Subroutine AXROT. A second modification mee 


volved the subroutine AXROT, which commuted the normalized 
direction cosines of the C-C bond axis using values of the 
direction cosines which relate the nuclei to the crvstal axis 
system as listed in section II. 

The vector, Cur representing that vector 
along the C-H, bond is expressed from values of the Hy tensor 


(Sectien fl) as 


x = Oey 
ieee OS 
Z, = Oe 47 eae 


Likewise the vector ent 1s represented by values of the H. 


tensor 
Xo = Oo 4c 
Va = SOO 
Z5 = -0.996 
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Addition and normalization of these vector comnonents yielded 
the direction cosines of the C-C bond with respect to the 
@avysetal axis system a*, b, c. 

(3) Subroutine AXR. As the methylene groun 
rotates through angle ¢', the hyverfine counling tensors, Hy 
and H, must change as a function of 9'. Hence, a further 
modification by the inclusion of subroutine AXR was required. 
The rotation of any vector a number of degrees ¢' about a 
vector (a,b,c) is represented by the transformation matrix 


resulting from evaluation of the following matrix multinlica- 


eLrons:: 
pb ac, \ yf hoes \ /b -a 
f A qa r cos¢ sing 0 Fi 7 0 \ 
P be | a Ale. = cum 
5 5 b | sind cos¢ 0 a | d 
t 


where d = hikes =a and an + b? + a = l. 
Assigning a, b, and c the values of the vec- 


tor components of the C-C bond axis as compvuted in AXROT, the 


resulting transformation matrix was applied to the H, and H 


1 g 
hyperfine tensors to yield the adjusted tensors as a function 
of ¢'. These values then became input to the basic SSESR 
program. 

(4) Subroutine REFH. As a final refinement, sub- 
EFOUtine REFH was written to insure orthonormalitv of tensor 
values, the result of which became inout to the AXR subroutine. 


With the modifications described it was pos- 


sible to calculate the stick spectra for any maanetic field 


oO 


Orientation and any desired angle of rotation, ¢'. From the 
Stick spectra the coupling constants for any angle, o', could 
be readily calculated by differences in the proper transition 
values. 

bam aeailat L Ont Ot Bist and <A? in 

Two additional modifications of SSESR were required 
for the evaluation of the average coupling constants as a func- 
tion of temperature. 

(1) Angular Devendence of Wave Functions. The 
solution of Schroedinger's equation for the harmonic oscil- 
lator as basically presented by Eisberaq [Ref. 5] alona with 
a Lagrangian formulation was used to find the system's eiqen- 
functions, but the displacement variable was $' rather than 
x as presented by Eisberg. 

The time indenendent Schroedinger wave equa- 


tion may be expressed as 


m2 3B gt) 
Strap OF a EV cas oe 


where V is the potential energy term, and I is the moment of 


inertia of the CH. fepate lye) - 


The potential energy may be expressed as 
leaner 
oe Ko e (A3-8) 


The kinetic enerdy as 


RN | SE he 2 et) i nN » 
ee xmr uw" = me (—) = x (> aa (A3-9) 


Therefore, the Lagrangian (LZ) for the system 
ic 


fee ne See (Ae =O) 


Solution of the Lagrangian yields the equation of motion 
g' - So" = 0 (A3-11) 
which upon solution yields the frequency, v, as 
2: les 7 
—_ xv = GAS es 


To determine the actual value of the fre- 


quency, v, the following Hamiltonian was studied, 





L2 
= g! Ki2s zs 
Be (A3-13) 
1 ee 
where = is derived directly from the classical mechanics 


expression for kinetic energy, 
KE = dmv 
2 
and L is angular momentum. 


Considering the angular momentum quantum 


mechanical operator, 





the Hamiltonian may be rewritten as 


_ © oe K,,2 
H — Sa Dore + xo (A3-14) 


which may be further rewritten using equation A3-12 as 
** 9.¢ 


ee 2 2 12 = 
H = sTpgre t 2m*Id , (A3-15) 


Operating on the wave function, w, with this 


Hamiltonian gives 


A? 3° PO") 2 2: 2 
which is an identity with equation A3-14. Hence, rather than 


determine the frequency, v, by such means as, sav the IR 


SD 


frequency between energy levels, it was possible to calculate 
K (and hence v) from equation A3-12 directly. 

The potential energy of the restricted rotor 
is expressed as 


V 


v= ma (l-cos 2¢') (A3-17) 


where V, is the barrier potential observed experimentally as 


5.8 kcal/mol. Expanding the cosine term as a series vields 
72 2) ee eo 
ay. Ve 7 7) 

+ higher terms. 
Neglecting higher terms, 

fe: 12 
V= V5¢ 

which is just the potential energy of the harmonic oscillator, 
so that 

= Prue, Lone = 

V = Vj¢'* = 5Ko (A 3=amee 


and 


sK = V, = 5.8 keete/mol,. (A3 =e 


In dropping higher terms and anproximating 
the harmonic oscillator the assumption has been made that 
Significant values of energy will occur only with a few de- 
grees either side of ¢' = 0 and that as $¢' approaches infin- 
ity, V approaches infinity as shown in Figure 6. 

Combining the results of equations A3-18 and 
A3-19, equation A3-7 is rewritten and transformed to dive 


3 Vege 
a + Ee Be gry = 0 - (3-20) 
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Harmonic 


Oscillator 
Potential 
Function 
and Wave 
Functions 
a 2 “a WA, - —— 
fo ~~" - = Vee ae 
Restricted Rotor oi ——__— 


Wave Function 


Figure 6. Comparison of restricted rotor and 
harmonic oscillator wave functions. 


Defining the terms 





ZIV. tie 
B= 6 (A3=2e) 
and ct = va ¢' (A3-23) 
equation A3-20 becomes 
3° 
ony) s 2 = 
Se a (B-a Z vegry = 0 


which with boundary conditions applied per Eisbergdg's proce- 


dures leads to the eigenfunction solutions 


Voc) = exp(-t°/2)H, (,) (A3-24) 
where An (c) are the Hermite polynomials for the resnective 
levels n = zero to n = infinity, n = zero representing ground 
state. It can be shown that these functions form an orthog- 


onal set over all space [Ref. 5]. 

Combining equation A3-21 with equation A3-23, 
and applying the proper units conversion factors, t was deter- 
mined to be 

C =Moewy OLS) (A3-25) 
where $¢' is in radians. 

Mites or oD abaglia ty Poe that a given state n 


will be found at an angle of CH. rotation, Oo 7 Ls 
P = rs 
n(c) ~ *n(z)*n (x) 


where V* (cys the complex conjugate of Yocy but since real 


wave functions are involved, 


and the integral of equation A3-3 takes the form 


co 2 
= - ei 2 ' 


= CO 


which when normalized, becomes 


90 
3.77629 -t* 9 F 
-90 


Limits of integration were chosen based on the assumotion 
previously mentioned that significant values occur only with- 
in a few degrees either side of ¢' = 0. This assumption was 
justified by considering that for the ground state where 


7 -¢* _ -ad'* 
* (ce ani 


After about 0.27 radians or 16 degrees, Pr has decaved to 


Zz) 
= of its original value. After 64° it is exvected that Ps) 
will have negligible values, hence integration over the lim- 
its of -90° to +90° was considered suitably close to those of 
toes PDlUS infinity. 

(2) Temperature Mependenice. sono cml OG celnel ras 
tion has been developed concerning temperature exvlicitlyv. 
The probability term derived above merely states the probabil- 
ity that a state n will be found at a given angle of CH. 
rotation ¢'. Boltzmann statistics, however, when applied, 
show the statistical distribution of electrons between the 
possible energy levels, the greatest percentage population 
occurring normally in the ground state level and decreasing 


exponentially with higher excited levels. 


og 


The population distribution in each energy 


level is expressed by the Boltzman weighting factor 


-E_/kT -E_/kT 


In& Si 


Ne = 7 aay = <j (A3-27) 


ge 
n 


where omy is the degeneracy factor, E is the difference be- 
tween zero energy and the fae state energy, k is the Boltz- 
man constant and Z is the partition function. Since all 
levels in the model are doubly degenerate, (Section V-B) dt, 
cancels in the equation. 

(3) Combining Coupling Constants, Angular and 


Temepeseinne Dependence (egueere mle <A> 


Equation A3-3 may now be rewritten 


-E_ /kKT 
( 90 ~r? - n | aa 
<A>, = see 3.77629A,, — Hy (cy (A3-28) 


Two subroutines were added to PREDICT II to 
perform the calculation. The first, AVGCC, evaluated the 
integral by numerical integration to yield the average cou- 
pling constants for each energy level, n, assuming 100% pvop- 
ulation of each level. The second subroutine, AVBZM, applied 
the statistical factor to yield the average coupling constants 


as a function of temperature, <A> These two subroutines 


AL 
appear in Appendix A as an integral vart of PREDICT II. 
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B. HYPERCONJUGATION EFFECT ON COUPLING CONSTANTS 

Hydrogen atoms attached to 8 carbon atoms (carbon atoms 
adjacent to the carbon atom having unnaired spin) undergo 
interaction with the unpaired spin. This interaction is usu- 
ally termed hyperconjugation and is a maximum when the carbon- 
hydrogen bond and electron orbital are coplanar. Hydrogen 
atoms attached to an a carbon atom with an unpaired electron 
in a hybrid orbital also undergo interaction due to orbital 
overlap. This can be loosely considered a tyne of hyper- 
conjugation and is the effect to be analyzed in this work. 

When the proton 2s orbital overlaps the carbon Dd, eb isa) 
there becomes a finite probability that the free electron 
will come in direct contact with the proton. This phenomenon 
is known as Fermi contact interaction and its effect is seen 
in the effective magnetic field term of equation A3-2. The 
Fermi contact interaction again changes the Hamiltonian val- 
ues and hence the coupling constants of the system. 

If the planar form of the CH,CO, radical is considered, 
a plane of symmetry exists and there will be no direct inter- 
action of the a carbon Pp, orbital and the hydrogen 2s orbital 
as seen in Figure 7a. Spin-spin is the sole interaction in 
this configuration. At any other angle we expect direct or- 
bital interactions as seen in Figure 7b, the amount of Fermi 
contact interaction being determined by the electron density 
on the protons. 

The input coupling constants used in SSESR were those for 


the planar model. In order to adjust these values for Fermi 


4) 


P, orbital 


“- C>Ce ae 


S eroital 


Figure 7a. CH. orbital configuration in the planar 


Z 


form of the CH,CO, radical 


Pp. orbital 





higuBpe: 7p. CH. Orbital Seonfiguration atter FOraci1on 
of the CH. group through angle ¢' showing 


Orbital overlap or mixing. 
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contact interaction (hyperconjugation), a molecular orbital 

Calculation [Refs. 4 and 19] using the CNDO approximation of 
Pople et.al [Ref. 16] was made to determine spin densities on 
the two protons for various angles of rotation, ¢'. As seen 
in Section V, a plot of spin density versus angle of rotation 


yielded a hyperconyjugation effect expressed by 


Ayy = Ay + ECO). Oo) seal SM, Deve) 
i 
where B = 1,420.4058 MHz, the hyperfine splitting constant of 
a hydrogen atom with a complete electron in the 2s orbital. 
This relationship was incorporated into the SSESR basic 
program as subroutine CCCA and used to adjust the computation 
of equation A3-3 for hyperconjugation as well as to determine 


the magnitude of the effect. CCCA appvnears as an integral 


part of PREDICT II in Appendix A. 


C. FLAP VIBRATIONAL MOTION 

An additional type of motion which can be termed as out- 
Of-plane flapping motion of the CH, group may be considered 
as a possible source of contribution to the temperature de- 
memaence of the coupling constants, (Figure 8). 

The effects of this type motion would be exvected to act 
very Similar to that discussed under the subject of hvper- 
conjugation except that in the flap motion, "bending" of the 
CH. group through an angle y is of interest. This action can 
be considered a further type of hyperconjugation effect lead- 
ing to mixing of the o carbon ee orbital with the hvdrogen 2s 


Orbitals but under a different mathematical relationship than 
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H 


Figure 8. Out-of-plane flapping motion of the 


CH, group in the CH,CO, radical. 


44 


that due to rotational motions As®with notational mortien, 

an increase of electron density on the protons (Fermi contact 

interaction) and a corresponding change in the coupling con- 

stants due to mixing of orbitals is expected for flap motion. 
Molecular orbital calculations based on the CNDO approx- 

imation (as used for rotational motion) were utilized to 

calculate the electron density for any angle of flan, y. 

From the results of the calculations, molecular energy val- 

ues were used to test the feasibility of this tvpe motion and 

electron spin density values to determine the functional de- 


pendence of the coupling constants upon the effect. 


D. SPIN RELAXATION EFFECT 
1. General Relaxation Theory 

To observe electron spin resonance absorption there 
must be a population difference between the two spin levels 
between which transition occurs. Carrington and McLachlan 
[Ref. 3] discuss this criteria in rather easy to understand 
terms. 

The rate of absorption of energy from the radiation 


field (power) is 


JE 
sz = N Pag (Eg-Ey) + NePay (Ey Eg) = sae (D1-1) 
where N is the population of eneray states E, OF Eg. Pip 


meenc probability of transition from state ao tors (ae) Same 
1S approximately equal to P. The difference in state povula- 


erem 1s n = Ny - Ne and follows a time dependent exponentially 


decaying function 


aes 
= nae 
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Equations D1l-1l and Dl-2 show that resonance absorption dis- 
appears at equilibrium. 

Smee NF = Ne in the absence of a maaqnetic field, 
application of a field must inevitably require interactions 
between nuclei and their surroundings as thermal equilibrium 
is approached. This interaction is called spin-lattice re- 
laxation and results from changing spin orientation as mag- 
netic energy is transferred to other degrees of freedom. 

In addition, other processes are known to exist in 
Spin systems due to nuclei-surroundings interaction which 
cause energy level shifts and hence an expected shift in cou- 
pling constants through solution of the effective spin Hamil- 
tonian of equation A3-l. These effects are called spin-spin 
relaxation, and through solutions of the Bloch equations 
yield correlation times representative of time the spin 
system remains in a given eneray configuration. 

2. Aveli Gatien of Reiexation Theory 

Through a qualitative analysis of the physical pic- 
ture of the rotating methylene groun, a relaxation effect was 
postulated. 

Consider as in Figure 9a the CH, methylene groun 
Oriented in a position having the maqnetic moments of protons 
Hy and H, precessing about effective magnetic fields renre- 
sented by vectors hee and ig consistent with the notation of 
Section A-3-a. When the two protons exchange sites as in 
Figure 96, each will see a different effective magnetic field. 


Tiere 1s then a finite probability that the wrotons will 
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oe Original Fields €¢. 
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Figure 9. Representation of the magnetic moments 


of nuclei Hy and H, fFlipp eng aor 


precession about one effective maqnetic 
field to another due to rotation of the 
CH. Group. 
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reorient or "flip" through an angle a into nrecession from 
about the original effective field to about the effective 
field now constituting their environment. The two limiting 
conditions occur when rotational frequency is very high corre- 
sponding to high temperature and when rotational frequency is 
very low - essentially nonexistant - corresponding to quench- 
ing of motion at low temperature. Alternatively, it may be 
said that at high temperature the frequency of rotation being 
fast will yield a high probability of flinv (rapid flip) be- 
tween the two states. At low temperatures, the probability 
of making a flip will be small, and approach zero at the 
point where rotation can be considered to be quenched. When 
the flip probability is low, significant relaxation effect is 
not expected. Hence a correspondingly greater average cou- 
pling constant is expected over that when flip probability 
pice elae Joy 

It should be noted that Figure 9 is only intended to 
present a representative picture of the relaxation effect. 
In a real system, the probability of flip or reorientation of 
the magnetic moment must be considered at every increment of 
CH. rotation. 

So. Catlemlaicron of Limiting "Conditiens 

In Section A-3-a Box's development for solution of the 
spin Hamiltonian was seen to yield transition frequencies and 
hence coupling Constants for the CH.CO5 radical. As derived 
from Box's relations the transition frequencies and hence 


coupling constant for low flip probability (low temperature 
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limit), are expressed by 


= | Raeeacael - | (D3-1) 


where the results of both fields are assumed oriented in the 
same direction at quenched rotation. Solution of this rela- 
tion was obtained by a slight modification to Crawford's 
[Ref. 4] SSESR computer program. 

Coupling constants were obtained for the ranid flinv 
(high temperature) limit through an approximation that at 
rapid flip conditions, the effective field experienced by 
each of the protons is an average of the two fields h, and 
h, experienced at quenched rotation. This approximation was 


expressed as 


— Thy ey) *P2 (4) | Tee 


which was solved, again, by the modification of SSESR to 
yield the coupling constants for the rapid flip limit. 

The difference between the limiting values vielded a 
relative shift in coupling constants for the relaxation ef- 
fect at any given orientation of external maqnetic field. 

4, Calculation of Intermediate Conditions 

Anderson [Ref. 1] has calculated frequency shift over 
a range of temperature but his calculations are based on an 
approximation that the motional Hamiltonian commutes with the 
time independent Hamiltonian, an approximation that does not 
apply to the model under consideration. To date a rigorous 
Calculation without using Anderson's or other approximations 


has not been developed (Freed and Fraenkel [Ref. 9]). 
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As a result, a modified solution of the Bloch equa- 
tions in the presence of a random field using Monte Carlo 
techniques was used to calculate the frequency and hence, 
coupling constant shift for intermediate values of flip prob- 
ability (temperature). 

A simple physical picture of the phenomenon calculated 
by the Monte Carlo technique is depicted in Figure 10 which 
represents a sphere with an orthogonal axis system X, Y, Z 
having its origin at the center of the sphere. At the origin 
is placed a proton. Vectors hes and Bey lie in the YZ pnlane at 
an angle a either side of the Z axis and represent the effec- 
tive magnetic fields hy and hs depicted in Figure 9. Assume, 
as an arbitrary starting point, that at the instant a maaqnetic 
field is turned on, the magnetic moment vector, M, of the pro- 
ton lies along the X axis with its tip falling on the sumeee= 
of the sphere at polar coordinates (85, do) - Let classical 
mechanics theory show that in the influence of a magnetic 
field, the vector M will precess about hy, tracing aspara 
along the sphere's surface from point one to point two. 
Finally, designate an increment of precession, A, such that 
after each increment a random decision as to whether or not 
to flip from precession about hy to precession about h, is 
made. The random decision is correlated to a selected pvrob- 
aod Lah Bee ion ce uslor, Po which upon problem solution becomes 
correlated with a value of temperature. As time progresses, 
the magnetic moment will first flip from precession about hy 


Le h, eve) ie back and forth as shown after increment two, four, 
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Figure 10. 


Pictorial representation of the trace of 
the magnetic moment of a proton as it 


flips from precession about effective 


magnetic field hy 1S hoe back and forth. 
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Six, Seven, vescght, ecteemmihesresult is a randomipath traced 
on the surface of the sphere. 

In the esr experiment an oscillating magnetic field 
along the X direction is observed. Hence, motion of the My 
component of the proton magnetic moment changing with time is 
the factor to be determined in the relaxation effect. If the 


probability of flip is zero, the M, cemponent will trace out 


x 
a cosine-like curve. However, when a flip occurs and M 


changes its direction of precession, the M,, component deviates 


X 
also from its original trace and begins tracing out a new 
cosine-like curve. As the M trace spreads over the surface 
of the sphere, it is expected that the My component will fol- 


low an approximate exponentially decaying cosine curve which 


can be represented by 


-a,t 
Ee) =e cosw.ast (D4-1) 


where ay is a decay constant, Wo is the frequency with zero 
probabilie, of flip, and ao is a frequency shift factor. 

Figure 1l is representative of the function, Fey: 
The Fourier transform of equation D4-1 (Figure llb) yields a 
Single peak of frequency w = W 9 a5 characteristic of the flip 
probability and hence temperature. Plots of flip vrobability 
versus w may then be made for a given orientation of magnetic 
field from which coupling constant shift may be correlated 
with temperature. 

a2) Monte Carlo Calculation 
To calculate this effect, a random (Monte Carlo) tech- 


nique of generating the most probable function, f'(t), 


Se 


M time 


Figure lla. A decaying cosine function 
representative of the X component 
of magnetic moment. 


ii ea (Power) 


w (frequency) 


Figure llb. The Fourier transform of Figure lla. 


approximating that of Figure lla was used. The comvnuter pro- 
gram was written by Tolles [Ref. 20] and appears in Appendix B. 

The magnetic moment M initially positioned at (8 7 Q) 
then flipped from about A 


was rotated initially about A to 


1 A 
rotation about ho, iack and forth aceording to, a first ordex 
Markofflan process. The problem is to find 
“ sg 
= M e 
Fw) Ss i - 

where Ty) is a function of frequency and the integral is a 
Fourier integral giving the spectral density. 

Let A be an angle through which the magnetic moment, 
M, rotates about ia prion to having the oppomtunity to fli 
from state h, to h.. The position of M after rotation by A 


il 2 


about h, is expressed by the rotation matrix maninulations: 


IL 0 0 cosA -sindA 0 1 0 0 sinscos¢ 

0 cosa sina Sink cosA ag 0 cosa -sina sinésing = 
' 0 -sina cosa 0 0 1 ae cahever = elefore) cos6 

cosh -sinAcosa sinAsina sin8cos¢ 

sinAcosa cosAcos*a+sin*a (-cosA+l)sinacosa sin@ésind 
-sinAsina (-cosA+1) sinacosa cosAsin‘*a+cos‘*a cosé 


which yield a new point (8, ,¢,) and hence a new value of the 
My component. 

Pecet scout lin, rotation of 4 about ii yields another 
point (6 7¢,) expressed as above with each sina replaced by 
-sina. The rotational transformation is made after each in- 


Sreneneaor FOtatron, A. 
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To perform the Monte Calo calewWilartmon a selected nun 
ber of points along the M, trace senarated by the value of A 
were chosen. For examvle, 36 points each separated by a A of 
10° would represent one cycle or one revolution of M about h, 
72 points would represent two cycles. A probabilitv, d, of 
having no spin flip after each increment A was chosen, a value 
of, for example d = 2 representing one chance out of two that 
eae flip will not occur after increment A (rapid flipping), a 
value of d = 100 representing Bo” chances “Gut on F100" Chatene 
flip will occur (slow flipping). A random number chosen com- 
pletely arbitrarily was fed to a random number generator 
(RANDU) which produced in turn another comoletely random nunm- 
ber representing a condition of either "yes fliv" or "no 
flip". RANDU was correlated to randomlv permit the system 
to flip in accordance with the desired flip probability 
Peor 100, etc. 

Sizar ting wach a “podnt Mo (Oy, do). rotation tarough 
A was calculated and a new value (65, >,) was determined. At 
A the random number generated by RANDU indicated the system 
Piped Or did not flip: If flip occurred, another value of 
(8, >) was calculated based on rotation about ie Lf nerear 
occurred, rotation was continued about ae After each incre- 
ment A, the decision was randomly made until the selected 
number of points was covered. Note that, especially for larce 
probability factors, such as 100, due to randomness of the 
generated yes/no decision, all points may be covered without 


any flip while on the other hand, several flins mav have 


Sy) 


occurred. To obtain a most probable picture of the true phys- 
ical phenomena, the calculation over the selected number of 
points was run over and over a selected number of times, the 
random number generator providing an ever increasing number 
of adjusted points along the MCE" (yy) curve. The final noint 
values were then fitted to the damved cosine function Fey? 
equation D4-l1, by a method of least squares. 

Output from the Monte Carlo calculation vields a trace 
of the damped cosine curve, the relative frequency, Wyss the 
exponential decay constant ay for the best fitting Eiey? and 


a listing of deviations of the calculated fr curve from 


t) 
the best fitting Fie) curve. 


eo  Sieloaloi Jinn? SE sey Ve Siem Cs Meneser etude 
At the time, t, when the first flipn decision is made 


(after A) the probability, Poe Gf wae Flip as 


i 
Lae 
After call decisions, one cycle has elapsed and time 
t = 2 seconds. 
Vv 
The probability, Poe of no flip has exvonentially 
decayed to 
po = pets (D6-1) 
n 
where 
es il 
Dee 1-4 : 
Hence equation D6-1 may be rewritten as 
-t/T . 
Pa =e (D6-2) 


where t is the time progressed and ive is a correlation time 
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representing the lifetime in a pamticulair ss eatemne Ole enema. 


ping to the alternate state. Using the relation 


ye @inD 
mime probability of no flip, at £ = = is 
ZT - i a i 
Rie exp ([— InD) = exp({ L/vt)) exp ( 1/t.). 
Hence, 
= ms 
"e ~ 2rinD nicer 


where a basis of frequency of 1 Hz has been used. 
The results of Tolles et.al. [Ref. 21] and Ohigashi 
and Kurita [Ref. 15] show 


a7 V/RT 


=| 
Fh 


(D6-5) 


where fo and V are the frequency factor and activation energy 
for the rotation of the CH. group. 


For given values of input parameters, A and d, as 
may be calculated by equation D6-4. From these values, tem- 
peratures are determined by equation D6-5 resulting in a cor- 
relation between flip probability, temverature, and frequency 
shift as calculated by the Monte Carlo computer program. 

ie HOumLeh RepresentatlonsOn theyepeerma 

The results of the Monte Carlo calculation yields 

valuable information about the predicted experimental spectrum 
representing the system. 


The symmetric Fourier integral of the curve, Fie)? 


best fitting the calculated points of f(t) is 


~ —-1iwt ~ay,t 
Fi = ‘ e e cosWyja,tdt 


Sy 


Solution of this integral yields, 


ae 
“ay a 
F “ieee (D7-1) 
1 O%2 1+=7 (w-wy a.) 


il 
The Lorentzian line shape function, characteristic of damped 


oscillatory motion and characteristic of esr absorption spec- 


tra, is 
oe) 
g = 2 a 2 (D7-2) 
(w) 1+T4 W=W 1 
where T, represents relaxation effects. Comparing equations 


D7-l1 and D7-2 the esr line shape parameters become 


Center frequency = Wy A5 


Ml 
to 
A) 


Half height = 2 
2 


f : —_— 
Line width = 7 = ay : 


These factors were used in analyses of the results of Sec- 


tion Y-F. 
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Iv. EXPERIMENTAL PROCEDURES 


A. CRYSTAL GROWTH 

As described by Crawford [Ref. 4] single crystals of zinc 
acetate dihydrate were grown by slow evaporation of aqueous 
solutions. A mass ratio of about 2.5:1l water to zinc acetate 
dihydrate was used to form a nearly saturated solution, ali- 
quots of which were allowed to evaporate from watch glasses. 


Single crystals formed in 12 to 16 hours. 


Bam erRYSTAL IRRADIATION 

Prior to irradiating a crystal, it was examined under a 
microscope to determine that the b axis (Figure 3) could be 
readily identified for proper orientation later in the esr 
Spectrometer. A single crystal was placed on scotch tape and 
mounted in the open window of an X-ray diffractometer. Three 
such crystals could be irradiated simultaneously. Exvnosure 
was made typically for about two hours at room temperature 
using 50 KV, 25 ma, and a copper target. Immediately after 
irradiation, the crystals were placed in closed containers 
and refrigerated until use. The acetate radicals formed by 
the X-irradiation remained stable under refrigeration of zero 
degrees Centigrade up to one month. The signal may remain 


even longer, but no irradiated crystals were kept longer than 


this before study. 
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C. CRYSTAL ALIGNMENT IN THE MAGNETIC FIELD 

In order to vary the orientation of the magnetic field, 
a crystal was rotated through angles 6 and ¢ as shown in 
Figure 12. The angles 8 and $¢ are converted to the direction 
cosines of the magnetic field with respect to a*, b, c axis 


system through the transformations: 


Axis Direction Cosine 
a* wees 0 

b -sin 8 sin $6 

Cc gin 6 cos > . 


Using these two angles it was then possible to describe 
any orientation in a three dimensional system. Since a 
Slight change in the magnetic field orientation will produce 
Significant change in the esr spectrum, precise orientation 
was necesSary. 

A crystal to be studied in the esr spectrometer was 
placed on a guartz rod having both ends ground flat and par- 
allel to each other, one end for crystal mounting, the other 
for clamping to an extension rod connected above the resonance 
cavity to an angle indicator for determining the angle 6. 
The crystal was held on the quartz rod by a spot of stopcock 
grease acting as a glue. This procedure held the crystal in 
place over the entire temperature range studied, -150°C to 
Tog Ce 

Prior to placing the crystal into the spectrometer, the 
angle ¢ was established. As mentioned, careful orientation 
was critical. To establish $6, a Bausch and Lomb microscope 


having a crosshair eyepiece over a calibrated revolving 
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Pigure 12... Single crystaleof zine g~cetate, dihydrate 
mounted on a quartz rod with the b axis 
coincident with the axis of the rod. 
Crystal is rotated by ¢$ and rod by 8 to 
establish magnetic field orientations. 
[Ref. 4]. 


platform was used. The quartz rod was placed on the revolv- 
ing platform and one crosshair was aligned pnarallel to the 
long axis of the rod. The platform (and rod) was rotated 
through the desired angle $¢ and the crystal was then aligned 
with the b axis parallel to the crosshair. The angle $¢ could 
be established to within +0.5° by this procedure. 

The angle 8 was established after the crystal and rod 
were placed into the spectrometer resonance cavity. The end 
of the rod opposite to the mounted crystal was clamned to the 
extension rod leading to an angle indicator positioned di- 
rectly above the cavity. This indicator established angles 


tO Wien ine O°, 


D. THE ESR SPECTROMETER AND TEMPERATURE CONTROL UNIT 

All data was taken using a Varian V-4502-13 esr spectrom- 
eter having a Fieldial V-Fr 2503 magnetic field regulator, a 
nine inch magnet, and a conventional reflection tvpe cavity. 
Temperature control was attained using a Varian variable tem- 
perature control unit which had been calibrated for proper 
NEEROgen slow to: wlthin ioc. 

During warmup and stabilization of the equipment, the 
angles ¢? and @ were established. With both angles set, the 
spectrometer was fine tuned at a magnetic field of (3210 + 
40.00) gauss, the 40.00 being the initial setting of the 
Fieldial. The desired temperature was established using 
liquid nitrogen and nitrogen gas per operating characteristics 
of the temperature control unit. Ten minutes or longer was 


allowed for the crystal and cavity to reach temperature 


equilibrium. 
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Most data runs were taken at a set orientation of 6 and 
d, annotated as (B) starting at a given temperature and vary- 
ing the temperature for each succeeding spectrum. Each time 
the temperature was changed, thermal equilibrium was allowed 
to be established prior to recording the spectrum. 

With all factors stabilized, first derivative spectra 
were recorded on an X-Y recorder over a 2.5 minute, 100 gauss 
sweep of magnetic field. Each absorption line of a spectrum 
was then determined and labeled using the field dial of the 
Spectrometer. To avoid backlash error of the field dial, all 
lines were approached from the low magnetic field side when 


identifications were made. Line values were read to 0.01 


gauss. Measurements could be repeated to within +0.03 gauss. 
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MV. RESULTS OF PROBLEM SOLUTION 


A. VERIFICATION OF SHIFT IN COUPLING CONSTANTS WITH 

TEMPERATURE 

The first objective in this investigation was to confirm 
the extraordinary shift in coupling constants with tempera- 
ture observed by Crawford [Ref. 4]. In Section IV, the meth- 
od of recording experimental data was discussed. Generally, 
data was taken at a fixed orientation of magnetic field (3) 
while the temperature was initially set at -150°C and then 
raised by increments of ten or twenty degrees, spectra being 
recorded at each temperature increment up to about 90°C. 
Above 90°C the crystal began to dehydrate and the esr signal 
due to the CH,CO, radical rapidly broadened, decreased in 
intensity, and finally disappeared, no signal being observed 
at 100°C. This effect was probably due to molecular reorien- 
tation upon dehydration leading to destruction of the CH,CO, 
radical. 

Figure 13 shows characteristic spectra obtained over the 
temperature range of -150° to +90°C. The particular spectra 
shown is that for a magnetic field orientation giving rise to 
two magnetically observable sites within the crystal. Two 
sets Of coupling constants are obtained from spectra of this 
type, one for orientation (3), the other for orientation (_g)- 
When $¢ 1s zero, the two magnetic sites are indistinguishable 


and spectra characteristic of that in Figure 1 are observed. 
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-150 °C 


70.5 °G 


FIGURE I!3 


Data as recorded from the esr spectra is tabulated in 
Table I. From the difference in values (in gauss) of the 
absorption lines, experimental coupling constants were ex- 
tracted for a corresponding temperature. 

Figure 14 shows characteristic plots of the coupling con- 
stants versus temperature for given orientations of magnetic 
field, while the total shift in coupling constants between 
low and high temperature are tabulated in Table II. 

For all orientations of magnetic field, the total observ- 
ed shift in coupling constants was in close agreement with 
that observed by Crawford, verifying that the large shift 
does in fact occur for the CH.CO, radical. Several addi- 
tional and very interesting phenomena were observed when the 
coupling constants were studied at increments of temperature 
aS in Figure 14. Crawford's cursory observation of the tem- 
perature effect was based solely on the coupling constants 
determined at two points, + 25°C and -140°C. This investiga- 
tion confirms Crawford's observation but in addition shows 
that for all orientations of magnetic field, the temperature 
effect 1s confined to a narrow range of about +30°C to -60°C 
as seen by the characteristic "S" shaped plots. A decrease 
in the temperature coefficient was to be expected at low tem- 
peratures as the system approached quenched rotation of the 
CH. group, however the decrease at high temperatures above 
+30°C was totally unexpected and eventually led to formulation 
of the relaxation model described in Section III. [Initially 


it was thought that the decrease in temperature coefficient 
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Figure 14a. Experimental results of the coupling 
constant dependence on temperature 
for selected magnetic field orientations, ( 
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Figure 14c. Experimental results of the coupling 
constant dependence on temperature 
for selected magnetic field orientations, 
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Figure 14d. Experimental results of the coupling 
constant dependence on temperature for 
selected magnetic field orientations, ( 
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mayaise duc 2@%a crystal legrapnicéhange in the system@ocenr— 
ring at about 30°C," However, "x Yay difiiracrien pattermncr or 
the crystal compared at approximately liguid nitrogen temnera- 
tures and at temperatures above 30°C showed that no such nhase 
change had occurred. 

The experimental data shows that for those magnetic field 
Orientations at which the two magnetically inequivalent sites 
can be distinguished there is a characteristic rise or drov 
in the coupling constants after leveling off as the temnera- 
ture is lowered. This drift of the counling constants remains 
unexplained. The total shift in coupling constants tabulated 
in Table II are based on the shift observed between the points 
where leveling occurs at low and high temveratures, the noted 
drift being disregarded. This procedure proved to be justi- 
fied based on calculations of the overall shift to be discus- 


sed under the relaxation effect results. 


B. MODEL DEVELOPMENT 

Experimental proof of the coupnling constant devendence 
upon temperature led to development of a model by which the 
effect could be explained. The work of Tolles et.al [Ref. 
21] and Ohigashi and Kurita [Ref. 15] showed rotation of the 
CH. group about the C-C bond; hence the restricted rotor was 
chosen as the initial model to be studied. 

The Hamiltonian for the restricted rotor has been discus- 


Sea in Section Lil “and was seem to be 


ae Z V = ' 
H= BLe =f 5 (1 cos2o ) 


vas 


TABLE Sil Sslotal coupling ecensmemteshifts obsexvedvexpem— 


mentally 
SD tteld Onitemeet lem Total Shift~ | 
| : (gauss) | 
pe 
| @.: : 0.35 : 
me) | 125 | 
| yt | 4.48 , 
| foe | 2.90 
Lo oe | 1.60 
3 
; : 
ey, | 1.86 
ae ae | 4.61 
| 45) | 
| \30/ | 1.80 
| — | 3.95 | 
: re. | 2a | 
: ( 60" | 4.28 | 
| arte | 2.50 | 
te | cae i 
ices) eae) | 
a 2.70 





Measured by the difference in extreme absorption lines. 


6 


where solution of the Hamiltonian matrix elements and comnut- 
erized diagonalization of the resulting secular determinant 
led to the eigenvalues characteristic of the restricted rotor 
system. Using the results of Tolles et.al. [Ref. 21] for an 
Bo@-H angle of llos9=vand™a C—-H bond lengewto: 1.08A, the 


mement of inertia, 2, of the CHA @sotation about seme we «spend 


2 
was calculated to be 1.69408A. The barrier potent Lal ay 

oe 5.8 kcal/mol was converted to units of V/H and foundere 

be 60.8524 THz (terahertz). Knowing I, B was calculated to 
be 0.298414 THz in units of B/h. From these values, the 
Hamiltonian matrix elements were readily calculated for 
integer values of m. 

As it was obviously impossible to work with a determinant 
of infinite size corresponding to the range of m values of 
minus to plus infinity, several determinants of increasingly 
larger size were solved until the variance in eigenvalues 
from one size to the next larger size was within 0.2% or 
aeeroximately .Ovl THz. It was found thatai29) x 29 deren— 
minant (m having values -14 to +14) was the smallest to meet 
these requirements. Tabulation of the eigenvalues for the 29 
mee? Geterminant appear in Table ZII1. 

Examination of the resulting eigenvalues for the restrict- 
ed rotor shown in Table III showed each energy level to be 
doubly degenerate as expected from the solution of Schroed- 
inger's equation which yields energy values 


mh? 


eT 





so that a -m or +m value yields the same eigenvalue. 
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TABLE ITI. Eigenvalues of the restricted rotor model: 






Energy Restricted Rotor 
Eigenvalues’ 
Tez) 






a ee 4.1853 
. A) ee 
i 
i 
j 


1 | | 2 eee 
Joes oe) teal 


mee? 20.2843 , 
: 7 20.2861 | 
| 


i 3 2s 221! 
| 27.8318 


4 | 34.9854 i 
35.0223 : 


41.7450 
41.8451 
48.2856 


aie is Ae F! 
54.0146 








j 
2 | 
| | 
| : 
i 
, 6 | 48.0342 


2 Te pm 


SS as 
Gore OZiag 


ee 


Values obtained using a barrier to internal rotation of 
5.3 Kealyme Ll. 
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Further examination showed that for the first four levels 
above ground state there was an avoproximately 8 THz difference 
between levels and that the ground state level itself was 
4.19 THZ2 whiclis approximately halt of ehemenewe  slevel di it— 
ferences. This observation suggested immediately that the 
restricted rotor system might well be approximated by an 


Hammonilc oscillator which has energy levels represented ie 


E = (n+5) hv (equation III-A2-1) 
mel WhLGh in itself ts one lime OF a reSEri cred tferor ele 


Harmonie oscillator approximation was applied to eehewrecrtr tee. 
ed rotor, and energy levels were calculated by using equation 
ITI-A2-1l as rearranged to obtain the eigenvalues in units 
consistent with those of the restricted rotor. (Equation 
ITI-A2-3, Page 27 ). Values of n equal Zero to eight were 
used in the calculation, the results of which were compared 
to the restricted rotor results and appear in Table IV. Ex- 
amination of Table IV shows that there 1s good agreement be- 
tween the two models in the first four levels above ground 
state, but that at higher energy levels, the oscillator model 
yields increasingly higher results. 

A Boltzman distributwon calculations teres =a0n appre amiate 
population of energy levels in the harmonic oscillator is 
shown in Table Vv. 

Based on negligible population of energy levels higher 
tian the fourth exerted State mw the harmonic Oservllacenr was 
assumed to be a good approximation to the restricted rotor 


and subsequent calculations were based on this approximation. 
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Harmonic Oscillator 
Eigenvalues 
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TABLE V. Energy level populations in the harmonic oscillator 





model. 

| | 

Energy Level Percentage Percentage j, 

| Population ; Population | 

. po=Aba2eanK a) Beek | 
| 

0 | 68.9 | 94.54 

1 : 23.5 | 5.15 

| | | 

| 2 | 5.8 | 0.29 : 

i ; 

3 | 16 | 0.02 i 
: ' 

4 Oy | 0.00043 | 


ene 
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Se oflFT OF COUPLING CONSTANTS DUE TO DIPOLAR. INTE RAGA HON 
In Section III development of the calculation used to 
determine the average coupling constant as a function of 


temperature (<A>,,) based on the effect of dipolar interaction 


his Giscussed. It was seen that the relation to be solved 
was 
@+90 -¢ mee 
<A> =: Soe H ————— do' ., 
soe he 3.77629A,, L - ae ‘ 


(equation III A-3-27) 
In the development of this equation, a number of comouterized 
subroutines were written to the basic computer vrogram SSESR 
to yield program PREDICT II which appears in Appendix A. Not 
considering hyperconjugation effects (Subroutine CCCA) which 
will be considered in the next section of this work, equation 
ITI-A3-27 was solved by PREDICT II to vield the shift in 


coupling constants for selected magnetic field orientations 
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Over a temperature range of +90°C to -150°C. The overall 
shifts are tabulated in Table VI. 

Before considering the significance of the results, the 
manner in which PREDICT II was utilized will be discussed. 

The theoretical development of each of the individual 
subroutines was discussed in Section III-A-3. Basic input 
data to the main program consisted of those parameters which 
Tolles et.al [Ref. 21] determined for the planar form of the 
CH.CO., radical plus parameters determined in this work to 
include, 

The number of protons in the radical = 2. 

The magnitude of the magnetic field or 

tie wuecrodt~ oem wEnequency =—13.6 THz: 

The principal values of the coupling tensors (coupling 

constants) as listed on page 18 of this work. 

The orientation of external magnetic field, (S) - 

The direction cosines of the H tensor as listed on 

page 20 %oneen1s wor: 

The first five eigenvalues of the harmonic oscillator 

as developed in section V-2 and appearing in Table IV. 

The temperatures at which the coupling constants were 

to be calculated. 

Thepilanits of~@imegraiaenjys—90°°toe+90ee 

The increments of d¢' used in the numerical integra- 

eno) a tee oar 


Miscellaneous constants as appropriate. 
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The output yielded the esr spectrum at any given angle of 
Fotation, $* paimeterms of traneaieemeerecUenetes sand 1 mten— 
sities as intermediate results. The coupling constants were 
then determined by frequency differences and automatically 
fed back into the program to yield the averaqe coupling seen 
stants <A>,, at given values of input temperatures for dinolar 
coupling effects. 

It was of prime concern to insure that PREDICT II was 
Pecating (Operating) properly and yielding weliablie) sesuiee. 
Hence, two operational tests were devised and used as checks. 

me. Lest One for Proper Rovamren 

The spectra computed at +90° rotation might be ex- 
pected to be coincident with that at -90°, and this would 
make an ideal check for proper rotation. Initially the re- 
sults of the two rotations did not agree which led to 
incorporation of subroutine REFH which insured orthonormality 
of the coupling tensors H, and H,. It must be noted, in 
addition, however, that the input coupling constants for pro- 
ton one are not exactly the same as those for proton two, 
hence the dipolar effect cannot be expected to yield iden- 
meal results .for +90° and -90° of sotation. To check for 
Ppoeeeer rotation, Hoth sets of anpueecouplingm@censtants were 
made identical (arbitrarily, identical to those of proton one) 
and rotation by PREDICT Li was found. to yieladeidentical seec— 
peowac +90° “anee-90° of rotation, “This prececume Cayemene 
good confirmation that PREDICT II was operating properly. 
For all future calculations except where indicated, input 
coupling constants were assigned their original values. 
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2. Test Two for Prover Rotation 
Using the transformation 


wNeSsn sine On 


= _ 
( Metniorsingecos? Ogg, M Ay 
oom, one 
an orientation of magnetic field oe) was rotated through 
an angle > of 90°, where the matrix M is that for the direc- 
tion cosines of magnetic field oe This transformation 


was performed with aid of a desk calculator and yielded new 
direction cosines for the magnetic field of 

Cronous 

roots? 

Ooo) TS 
These direction cosines were then fed to PREDICT II as input 
direction cosines for the magnetic field while also setting 
the input coupling constants identical as in test one, and 


the spectra were computed for angles of CH, rotation of -90° 


2 


and +90°. If the program was rotating proverly, the resulise 


ing spectra were expected to correspond to that spectrum for 


180 
0 


of the CH, group. Agreement was excellent and it was conefud- 


an input of magnetic field orientation ( ) at O° rotati@nm 
ed that PREDICT II was rotating properly. 

Table VI shows that a shift in coupling constants be- 
tween -150°C and +90°C did occur due to dipolar coupling, 
however the magnitude of the effect in no instance was that 
ObServyer=arperamencally. The direction of the dipolar cou- 


pling effect was seen to be dependent on the external magnetic 
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TABLE VI. Calculated coupling constant shift over tempera- 
ture rahge of -150°2¢ to. 90" 6 sdiew tendo. 
coupling.’ 






Riedie. Orilentaivon Shift ae 






ae 0.013 

: “ay | +0.006 : 
| (60 0.000 ! 
4 +0.033 : 
22) -0.035 : 
| 3 : -0.054 : 
: 45) +0.045 
| i : =O. Oo : 
| ec i = (0) sone . 
: | 
| a +0.048 

__<, | 40.023 

eu 0.013 
| ew -0.079 | 
| oa } -0.119 | 


All shifts are determined by subtracting the coupling con- 
Beant ato lo0°%Cetromuthiatwat +9026- 
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Pied laos ic hovcwthatescomeemagnetic fields caused an 
effect which led to an increase of coupling constants with 
increasing temperature (dipolar coupling values positive) 
while others led to a decrease in coupling constant with in- 
creasing temperature (dipolar coupling values negative), the 
general effect observed experimentally. This effect was also 
observed graphically with plots as shown in Figure 15 where 
the calculated coupling constants due to dipolar interaction 
have been plotted against angle of CH. rotation, do tess 
several magnetic field orientations. 

As the temperature is raised, greater angles of @' 
are attained. Hence, analysis of Figure 15 showed that in- 
creasing temperature would be expected to cause an increase 
in coupling constants for those plots having positive curva- 
ture at @' equal zero (e.g. the Gres) plot) and a decrease 
in coupling constants for those plots with negative curvature 
at ¢' equal zero (e.g. the (a plot). The curvature Of@ea. 
plots and the direction of coupling constant shift in Table 


aa Orientation showed curvature of nearly 


VI agreed. The ( 
zero and this too agreed with the nearly zero value seen in 
Mable: VI: 

These results showed that although some shift in cou- 
pling COnseantewoccunred due to dipolar interaction, the 
direction of shift was not always consistent with that ob- 


served experimentally and hence dipolar interaction was not 


the primary cause of the effect. 
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It was of note that although the calculated shift was 
small, the magnitude was the same order as that observed for 


the methyl radical by Zlochower et.al [Ref. 23]. 


De SHLET OF COUPLING CONSTANTS DUE TO FERMI CONTACT INTERES 

ACTION 

The theory of hyperconjugation and resulting Fermi contact 
interaction was discussed in Section III-B. 

Using a computerized orbital calculation based on the CNDO 
approximagion [Ref. 19] the proton 2S molecular orbital coc 
ficients were calculated for various angles of rotation, $'. 
The coefficients squared yielded the amount of free electron 
Spin density that was imparted to the proton 2S orbitals for 
any angle of rotation. The solid curve of Fiqure 16 shows a 
plot of the resulting spin densities versus the angle of rota- 
tion, >', which upon examination appeared to closely approx- 
imate a function represented by 

Ean) = Asin*29' 
Of amplitude gAg= 0.273 x oo The dotted curve of Figure 16 
zs a plot of @his function to show that, in fact, the@aeipen= 
conjugation effect did very closely approximate the sin‘*29' 
relationship. 

If the entire electron (spin density = 1.0) were found Gam 
apLOecOn 25 Orbeecall, then the coupling constants of the plamean 
radical (where due to orbital symmetry, no spin density an- 
Peds see prOtOmeesmorbital) would be increased by «the 
factor B = 1,420.4058 MHz, the hyperfine coupling constant of 


PeIeEOden atom with a complete electron in the 2s orbital. 
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Hence, the input coupling constants for the planar radical, 
A,, were adjusted for Fermi contact interaction through the 
relation 
= 2 oe 
aa = A. + B(0.273x10 “sin” 26') 
This relation was incorporated into PREDICT II through sub- 
routine CCCA (Appendix A) and values of the average coupling 


constants <A> were again calculated having both dipolar and 


-! 
Fermi contact interaction incorporated. Column 2 of TablG@iame 
shows the results obtained for the total shift in coupling 
constants over temperature range -150°C to +90°C which can be 
compared with column 3, the values for dipolar interaction 
only. 

Since the hyperconjugation effect has no dependence on 
Orientation of magnetic field but only on angle of rotation, 
o', it was expected that the effect should be constant for 
all field orientations. By subtracting the dipolar effect 
from the total effect, the magnitude of hyperconjugation 
effect on the coupling constant shift was determined, the re- 
sults being tabulated in column 4 of Table VII. As expected, 
the effect 1s essentially constant for all magnetic field 
orientations having a value equal to about -0.19 gauss over a 
temperature range of -150°C to +90°C. The negative value 
corresponds to a decrease in coupling constant with increas- 
ing temperature. 

The cabveulated average coupling constants incorporating 
both effects were plotted versus temoerature giving plots 


SieeaGteristic of Figure 17 for all magnetic field orientations 
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TABLE VII. Calculated hyperconjugation effect (Fermi con- 
tact interaction) on coupling constant shift. ° 





| Shi fsewith Shite waem Shift Duemees 
Field Dipolar and Dipolar Hypercon)jugation 
lorientation Fermi Contact Interaction (Fermi Contact) 
| | Interaction Only | 
| | (gauss) (gauss) (gauss) 
! | 
| ve 0.009 =O Sa =()LOE 
30 | 
: 45 : : J | 
) 60 : Ore 0.006 i O Jha | 
i @ 10,196 0.00 i aout ! 
60 e e e | 
} ; x i { 
; {75\ . | . - ! 
i \30/) | OF las HO OOS : OAKS ul, ! 
! . | : 
) ee t =-0.220 -~0.035 : -0.185 | 
; 60 t j ‘ } 
| 90° ! : | 
| a OS -0.054 | =i().5. 1.9 1 | 
: : | 7 
—— EE E>E————— 
Sy | | ! 
! 90 \ | | 
7 75) -0.283 -0.091 | -0.192 : 
90. 
! 90/ = OleZ Ore, -0.079 : -0.185 ; 
{ 
{ 
| 27) eZ o +0.048 =O eli 
x oO } f 
| 135 : | 
.) -0.167 +0.013 , ~@% 180 | 
. { 
: ‘ | 
1 = =6-0.013 ) -0.185 | 
| 
| | 
(165) De2G8 -0.079 | -0.189 3 
| 
C8°) -0.311 +0.119 | -0.192 
i 





i shifts are determined by subtracting the counling 


cCconstantwea tua! 50° Cammaommtthate at +90°C. 


S)all 


studied. Examination of these plots and the data of Table 
VII led to several conclusions regarding the model as calcu- 
Kaeed ~O this point. 

First there appeared to be a definite effect on the coue 
pling constants when both dipolar and Fermi contact interac- 
tion were considered. The overall effect led to an increase 
in coupling constants with decrease in temperature as observed 
experimentally, the primary effect arising from Fermi contact 
interaction. The plots of Figure 17 indicate a decrease in 
the temperature coefficient as lower temperatures were av- 
proached indicative of quenching of the rotational motion. 
This observation agreed with experiment but it was of note 
that no such temperature coefficient decrease occurre@iasr 
higher temperatures in the calculated effects. Of greatest 
Significance, however, was the magnitude of the calculated 
effects which was seen to be only about 0.1 to 0.3 qauss, 
not nearly that observed experimentally (1-6 gauss). 

Figure 18 compares a typical experimental result with that 
calculated on the basis of dipolar and Fermi contact 
interaction. 

Based upon the foregoing observations, it was concluded 
that although the calculated model showed an overall effect 
generally of the form observed experimentally, and although 
dipolar and Fermi contact interaction can be expected to con- 
tribute to the shift in coupling@eonstants, these were Goma 
a small part of some larger effect which was yet to be 


established. 
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Figure 17. Dependence of coupling constants on tempera-— 
ture due to dipolar and hyperconjugation 
effects calculated. 
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Calculated dependence of coupling constants 
on temperature due to dipolar and hyper- 
conjugation effects compared to experimen- 
tally observed temperature dependence. 
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E. SHIFT OF COUPLING CONSTANTS DUE TO FLAP MOTION 

Out-of-plane flapping motion was discussed in Section 
III-C and was seen to be a type of hyperconjugation effect 
leading to Fermi contact interaction as with the rotational 
model but under a different mathematical relationship. An 
attempt was made to determine if this type vibration was 
Secucrinagwim the CH, CO. radical and the magnitude thereof. 

A molecular orbital calculation of the same type used 
for the rotational model [Ref. 19] but with various angles 
of flap led to the proton 2S orbital spin densities tabula- 
ted in Table VIII. In addition, the relative energies of 
the radical at the various angles of flap were also deter- 
mined from the molecular orbital calculation and appear in 
Table VIII. 

No simple function relating the coupling constants to 
angle of flap could be obtained from the data of Table VIII. 
However, a plot of the energy versus angle of flap as seen 
in Figure 19 showed a double minimum which indicated that in- 
deed, a flapping vibrational type of motion may exist within 
ene system. The fact that the emergy offthe radical de- 
SmeaseS GO .a Minimum at about S205 of flap vindicated that if 
Flap does exist, then the most stable form of the radical may 
not be planar as reported by other workers. [Refs. 15 and 21}. 

A means by which this possibility could be tested would 
be to lower the temperature of the crystal to liquid helium 
values in an attempt to freeze out the two stable forms, one 
with an angle of flap of about + 20° the other of - 20°. An 


esr spectrum would then be expected to disclose four 


aS 


Tig Vila results of molecular orbital calculations fox 


flap vibrational motion. 


Angle of Pere Ho. Spin Density” Relative Energy 
(degrees) x 10 of Radical (aeu) 

OF 000 =o 2] oC 
0.004 -52.10980 
0.049 =52,L0go0 
O02 0 =5 2%, LUG 
0.182 =52. 11056 
ae Jesh) -52.11097 
OC) pam ie -52.11110 
OG ake =S2eL Pte 
02096 =o 2. L Oe 
Oo -52.11081 
0.054 -52.11059 
0.040 -52.11030 
00:30 -52.10994 
0. 00% =52. 107% 
O(a e) I =52,.0 50m 
Or 002 -32.095% 
Ore 0'0 =52.0000 
Ore2 0 -52.07244 
O21 29 =52. 05000 





IU ; : 
Negative angles of flap yield nearly identical results. 


2 ; ioe : 
Values reported to three significant figures. 
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inequivalent sites rather than the two observed at the temper- 
atures studied in this work. If the four sites were to be 
identified, then raising the temperature through increments 
would yield data from which, via Arrhenius plots, a barrier 

to internal flapping motion could be determined. 

Flap vibration is not the most likely mechanism to ex- 
pect. No experimental evidence of the effect was observed 
over the temperature range of -150°C to +100°C. Based on the 
more obvious rotational mechanism which can be substantiated 
by experimental data, it was concluded that if flapping does 
occur, it is most probably characterized by a very low fre- 
quency and large amplitude with insignificant effect upon 


the system's coupling constants. 


F. SHIFT OF COUPLING CONSTANTS DUE TO RELAXATION EFFECTS 

The development of the relaxation model and methods of 
calculation were discussed in Section III-D. It should be 
noted that the development was based again on the original 
model of the CH, group rotating about the C-C bond of the 
CH.,CO, radical. 

The first step in making the calculation was determina- 
tion of the limiting conditions and hence the overall shift 
in coupling constants between these conditions. SSESR with 
evo lligittemedrrication was Used to calculate the Vimits- 

For a given magnetic field orientation, the relative 
coupling constants at the low temperature limit (low flip 
probability) were obtained through solution of Box's [Ref. 2] 


relation, 
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+ hy | 


the results of which appear in Table IX column 2. 

The high temperature limit was based on the anproximation 
that at high temperature conditions, the rate of flip is so 
rapid that the effective field experienced by each proton will 
be an average of the two fields hy and hs as expressed by the 


relation 


hy (_) | 


a 
Results from solution of this relation are seen in column 3 
of Table IX. 

The difference between the two limiting coupling constants 
yielded the overall shift 4. (calculated shift). VSnncewthe 
program calculated the values in units of megahertz it was 
necessary to divide the resulting overall shift by the con- 
version factor 

2.80 MHZ/gauss 
in order to compare the results with experimental overall 
shifts (A). Column 4 of Table IX shows the final calculated 
shift and column 5 the overall shift observed experimentally. 

The values of A were plotted versus the values of Aa 
Semcecl in Figure 20, the points clmster about apcreeaicle 
line, deviations being assessed to other effects which have 
not been isolated. 

Each point on Figure 20 represents a different maqnetic 
field as tabulated in Table IX. It is of note that approx- 
imately ten hours of experimental work and computations were 


Fequi bed, EOue@ctablici! One point. It was moted,in Sectrenm /-4 


Og 


TABLE IX. 


Comparison of overall calculated coupling constant 
shift with experimentally observed shift. 
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that experimental plots of coupling constant versus tempera- 


6 
4) 


coupling constants with lower temperatures after they leveled 


ture for those cases where $¢ # 0 (°%), showed drift in the 

emt at about -GO°C. The valuesmof A. used in Table IX and 
Figure 20 were, for these cases, measured between the high 
and low leveling off points, and as Seen in Figure 20 agree 
well with the calculated shifts. Initially, the deviations 
from linearity of Figure 20 were thought perhaps to be due to 
the drift observed experimentally. However, examination of 
the data showed no general correlation to confirm this fact. 

Based on the linearity of Figure 20 it was felt that the 
method of calculation for the limiting conditions was repre- 
sentative of the system. Having determined the limiting 
conditions and hence the overall shift in coupling constants, 
tae calculation of coupling constants) for intermediate flip 
probability (intermediate temperature) was undertaken using 
the Monte Carlo technique discussed in Section III-D. 

The Monte Carlo technique is based upon random statistics, 
hence the more points on a given curve that are calculated 
and the more times each point is calculated, the more reliable 
the final results are expected to be. Numerous calculations 
of the damped cosine function (Equation III-D4-1, Page 52 ) 
approximating the M, component of magnetic moment and the 
resulting frequency were made by varying the input parameters 
Of flip probability, d, (from 2 to 100), precession increment, 
A, (from 2 to 10), number of cycles calculated (1 or 2), ran- 
dom number input, and the number of times a particular set of 


parameters was recalculated (100 to 400 times). The angle a 
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Figure 20. Comparison of total coupling constant shift 
calculated (A) with that observed exveri- 


mentally (A). 
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between the Z axis and the effective magnetic field was chosen 
as 25.84° which permitted sfrequene Bosse (nepmesented acm ie eae 
tive values between 0.90 and 1.00 for a given flip probabil- 
ity. The magnitude of the results were independent of a, but 
the specific value was chosen for ease in converting the fre- 
quency data to increments of overall shift in the coupling 
Constants. 

For each value of d and A used to determine the best fit- 
ting damped cosine function, values of correlation time, Tae 
were calculated by equation III-D6-4 based on a frequency of 
precession, v, Of one hertz. The values of oa calculated 
thusly were plotted against the frequency of the best fitting 
damped cosine curve with the results shown in Figure 2l. 

There is a considerable scattering of points in Figure 21 

due to the random statistical method of calculation. However, 
it was not difficult to describe a curve most representative 
Of the greatest majority of points. The "S" shaped curve 

was noted to be most characteristic of that type function ob- 
served experimentally. 

For further calculations, only those points lying on or 
close to the curve described in Figure 21 were used, these 
conditions being considered as most representative of the re- 
laxation model. Examination of the spectra calculated by 
PREDICT II for given temperatures and magnetic field orienta- 
tion showed that the average frequency of precession exveri- 


enced by the system was about 30.28 MHz. Hence, each 
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Figure 21. 
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Relative frequency as a function of correlation time based on 
correlation time calculated on a1 Hz basis. 
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correlation time as calculated based on one hertz was adjusted 
Ome cer consistent with the real system where 
= OD ole. 
‘real ie 


Using equation III-D6-5 (Page 57 ) and 1 values of 


real’ 
temperature were calculated to correspond to the frequency 
for a given probability d, the results of which are tabulated 
mao Table xX. 

The plot of relative frequency versus temperature in Fig- 
Mire 22 shows that the "S" shaped function Gharacteristic of 
the observed experimental function was retained. 

The final step was to convert the data of Figure 22 or 
Table X to coupling constant shift for specific magnetic 
field orientations. AS Seen in Figure 22 the frequency was 
represented as a relative frequency shift between 0.90 and 
1.00 where at any given temperature the percentage of total 
shift between the limiting frequencies and hence counvling 
constants could be readily determined. For example, ata 
temperature of 0°C, 43% (0.943) of the total shift in cou- 
pling constants is expected to have taken place. Since very 
good agreement was obtained between the experimental and 
calculated low temperature limits and overall shift of cou- 
pling constants, the percentage shift at a given temperature 
was applied to observed overall shifts and curves were plot- 
ted based on the calculations just discussed, the results 
being compared with experimental curves as seen in Figures 
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TABLE X. Results of Monte Carlo calculations correlated 
tLOeLeCMperaLuLre. 


(Probability Relative Frequency Temperature (°C) 
[Factor (d) Sais t 

2 0.9010 +64 : 
| z 0.9007 | +93 ! 
! 3 0.9021 +44 | 
5 | 0.9107 +25 | 
| 6 | 0.9195 ! +19 | 

8 ! 0.9327 : +10 . 
10 | 039422 | + 4 : 
| oats | 0.9600 | = 4 | 

18 | 0.9725 | -12 7 

25 | 0.9840 ake : 

30 : 0.9863 | -23 

32 : 0.9352 | + 9 ; 
| 50 | 0.9794 | ik . 
| 60 i 0.9915 37 | 
75 0.9899 ye | 


100 0.9914 -34 | 
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Figure 22. Relative frequency as a function of 


temperature based on an average 
precession frequency of 30.28 MHz. 


ILIGhY: 


43 


42 
ww 
3 
fe 4l 
a 
Ae) 
G 
10) 
ny 
un) 
SG 
5 
43 
m 
G 
4 
_ 
OQ, 
= | 
O 
OU 
42 
4l 


Eo 


Figure 23a. 
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Comparison of the calculated and exnerimental curves was 
made by noting the temperatures at which the temperature coef- 
ficients changed sharply (temperatures where the curves level 
off) and the inflection points. 

The calculated curves followed a function closely approx- 
imating the experimental curves, but with a shift upward in 
temperature of 10°C to 20°C. The position of the calculated 


curves depended heavily on the value of T shifting down- 


real’ 
ward in temperature for lower values of average frequency 
experienced by the system. However, this fact was of minimal 
value in that from calculated spectra results, the average 
frequency of 30.28 MHz was the proper value to use for the 
system. 

Although the calculated curves did not agree with experi- 
mental curves as well as would have been desired, based on 
the random statistical method used in making the calculation, 
it was felt that the results obtained were in fact very good 
in that it was shown that the relation observed could be 
essentially constructed entirely from theoretical calculations. 

Of greatest significance, however, was the fact that the 
temperature dependence observed in the GHAGe, radical esr 


spectrum could now be explained as being due primarily to a 


relaxation effect as was developed in this work. 
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Vibes CONCLUS LONS 


As originally noted by Crawford [Ref. 4] a marked temper- 
ature effect upon the esr spectrum of the Seo radiedt 
formed in X-irradiated zinc acetate dihydrate was verified. 
Decrease in temperature corresponded to increase of coupling 
constants extracted directly from the spectra of a given maaq- 
netic field orientation. The total shift in coupling con- 
stants does not spread over the entire temperature range of 
-150° to +90° studied but was observed to occur between the 
rather limited range of -60°C to +30°C. 

Development of a mathematical model to describe the ob- 
served effect led to considerations of several phenomena 
expected to take place in molecular systems, the results of 
which indicated that the effect was due to a combination of 
factenes 

Dip@far interaction due to rotation Mf the CH. group about 
the C-C bond agrees closely with that observed by other work- 
ers £0605 a rotating methyl group but does net account folgene 
observed experimental effect. Fermi contact interaction due 
to hyperconjugation effects in the rotational model shows cou- 
pling constant shifts with temperature of the general nature 
observed experimentally, but again, the magnitude of the effect 
accounts for only a very small portion of the observed shift. 

Flap vibrational motion, although perhaps taking place in 


the system as indicated by molecular orbital calculations, is 


a2 


not observed experimentally and is not expected to be a fac- 
tor of any great importance on the coupling constants. 

The primary effect causing the observed temperature denpen- 
dence is due to a relaxation mechanism within the molecular 
system. This effect is confirmed experimentally and a model 
has been developed using random statistical calculations. 

It is felt at this time, that given an orientation of exter- 
nal magnetic field, and a temperature at which the CH,CO, 
radical is stable, the coupling constant can be calculated 
fPemwithin t1.0 gauss. 

The procedures used to studv the temperature effect are 
not restricted solely to the Zinc acetate system but could be 
applied to other systems in which similar relaxation proc- 
esses are expected to occur. 

Although it appears that this work combined with the work 
of Tolles et.al [Ref. 21] and Ohigashi and Kurita [Ref. 15] 
nearly exhausts the CH,CO.,, radical system, there are yetre 
few remaining studies which would be worthwhile to undertake. 
Flapping motion has not been confirmed nor proved to be non- 
existant. The drift in coupling constants at low temperatures 
in some orientations of magnetic field has not been explained. 
Each of these two factors deserves an explanation but of most 
value to this work as well as to a general contribution to 
the field of molecular studies would be the developvment of a 


rigorous mathematical model by which the coupling constants 


could be calculated over the entire range of temperature. 


eS 


This is a challenge yet to be met in the field and it is the 
hope of this author, that this work may be a major step in 


that direction. 
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APPENDIX A 


A computer program, PREDICT II, for calculating taegay- 
erage coupling constants of the CH,CO, radital as a function 
of temperature. Various combinations of the subroutines may 


be used to calculate intermediate results Mas was) discussed 


in the foregoing thesis. 
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SUBROUTINE AVBZM(AAVGT1AAVGT2,EN,T, FAVE) 
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SUBROUTINE CCCA(AF,ANGLE,A,AA) 
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APPENDIX B 


A computer program to calculate spin relaxation effects 


Hsing a Monte “Carlo ™=reciiuque. 
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